Ongoing declines for the world's amphibians in the face of emerging threats.

Systematic assessments of species extinction risk at regular intervals are necessary for informing conservation action1,2. Ongoing developments in taxonomy, threatening processes and research further underscore the need for reassessment3,4. Here we report the findings of the second Global Amphibian Assessment, evaluating 8,011 species for the International Union for Conservation of Nature Red List of Threatened Species. We find that amphibians are the most threatened vertebrate class (40.7% of species are globally threatened). The updated Red List Index shows that the status of amphibians is deteriorating globally, particularly for salamanders and in the Neotropics. Disease and habitat loss drove 91% of status deteriorations between 1980 and 2004. Ongoing and projected climate change effects are now of increasing concern, driving 39% of status deteriorations since 2004, followed by habitat loss (37%). Although signs of species recoveries incentivize immediate conservation action, scaled-up investment is urgently needed to reverse the current trends.

Middle Jurassic fossils document an early stage in salamander evolution.

Salamanders are an important group of living amphibians and model organisms for understanding locomotion, development, regeneration, feeding, and toxicity in tetrapods. However, their origin and early radiation remain poorly understood, with early fossil stem-salamanders so far represented by larval or incompletely known taxa. This poor record also limits understanding of the origin of Lissamphibia (i.e., frogs, salamanders, and caecilians). We report fossils from the Middle Jurassic of Scotland representing almost the entire skeleton of the enigmatic stem-salamander Marmorerpeton. We use computed tomography to visualize high-resolution three-dimensional anatomy, describing morphologies that were poorly characterized in early salamanders, including the braincase, scapulocoracoid, and lower jaw. We use these data in the context of a phylogenetic analysis intended to resolve the relationships of early and stem-salamanders, including representation of important outgroups alongside data from high-resolution imaging of extant species. Marmorerpeton is united with Karaurus, Kokartus, and others from the Middle Jurassic-Lower Cretaceous of Asia, providing evidence for an early radiation of robustly built neotenous stem-salamanders. These taxa display morphological specializations similar to the extant cryptobranchid "giant" salamanders. Our analysis also demonstrates stem-group affinities for a larger sample of Jurassic species than previously recognized, highlighting an unappreciated diversity of stem-salamanders and cautioning against the use of single species (e.g., Karaurus) as exemplars for stem-salamander anatomy. These phylogenetic findings, combined with knowledge of the near-complete skeletal anatomy of Mamorerpeton, advance our understanding of evolutionary changes on the salamander stem-lineage and provide important data on early salamanders and the origins of Batrachia and Lissamphibia.

A new species of Bolitoglossa (Caudata: Plethodontidae) of the Bolitoglossa franklini group from an isolated cloud forest in northern Guatemala.

We describe Bolitoglossa qeqom sp. nov. from an isolated cloud forest in Cerro Guachmalén, Alta Verapaz, Guatemala, based on multiple lines of evidence (morphological, molecular, and biogeographic data). This region comprises a mountain ridge without previous herpetological surveys. The new species is a large salamander with uniform purplish-black coloration and is distinguished by having relatively long legs with only one costal groove between adpressed limbs, numerous maxillary teeth, few vomerine teeth, only one phalange free of webbing in digit III of feet, and a relatively short tail. It is geographically closest to its sister clade of B. lincolni + B. franklini and the xeric Chixoy river canyon appears to be the major biogeographic barrier that isolated the new taxon. The cloud forest inhabited by this species has undergone severe habitat destruction in the region and land conservation actions are urgent.

Taxonomic clarification and neotype designation of two Taiwanese salamanders (Amphibia, Urodela, Hynobiidae).

Taxonomic identification of two Taiwanese salamanders, Hynobius sonani (Maki, 1922), and H. formosanus Maki, 1922, was clarified. Based on the descriptions and measurement data in the original description, we confirmed that H. sonani has been treated as H. formosanus and vice versa, which might have been caused by confusion of sampling localities of the two species. Because we could not find the holotypes of the two species and concluded these were lost, we herein assign neotypes for these species for precluding future taxonomic confusion.

Amphibians in Zootaxa: 20 years documenting the global diversity of frogs, salamanders, and caecilians.

Zootaxa is a mega-journal that since its inception, 20 years ago, has contributed to the documentation of the planet's biodiversity. Its role concerning terrestrial vertebrates has been crucial especially for amphibians, which are the most threatened class of vertebrates. As current editors of the Amphibia section, we reviewed the state of knowledge of taxonomic publications on amphibians over the last two decades (from 2001 to 2020). Our review reveals that 2,533 frogs, 259 salamanders, and 55 caecilians have been named in these 20 years, mainly in the tropical regions of South America, Asia, and Africa. More than half (57%) of these species descriptions were published in only 10 journals. At least 827 species of the new amphibians (29% of the total) were described in Zootaxa. This mega-journal has served also as a place of publication for monographs and systematic reviews, in addition to short articles documenting the vocalizations of anurans and the morphology of embryos and larvae. Its efficient evaluation process, the freedom of manuscript length, including full-color figures, and free of cost for the authors, has made Zootaxa a favorite for amphibian researchers. In an era of accelerating rates of biodiversity loss, documenting, describing, naming, and proposing evolutionary scenarios for species is, more than ever, an urgent task.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra.

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders.

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Comparative osteology of the hynobiid complex Liua-Protohynobius-Pseudohynobius (Amphibia, Urodela): â . Cranial anatomy of Pseudohynobius.

Hynobiidae are a clade of salamanders that diverged early within the crown radiation and that retain a considerable number of features plesiomorphic for the group. Their evolutionary history is informed by a fossil record that extends to the Middle Jurassic Bathonian time. Our understanding of the evolution within the total group of Hynobiidae has benefited considerably from recent discoveries of stem hynobiids but is constrained by inadequate anatomical knowledge of some extant forms. Pseudohynobius is a derived hynobiid clade consisting of five to seven extant species living endemic to southwestern China. Although this clade has been recognized for over 37 years, osteological details of these extant hynobiids remain elusive, which undoubtedly has contributed to taxonomic controversies over the hynobiid complex Liua-Protohynobius-Pseudohynobius. Here we provide a bone-by-bone study of the cranium in the five extant species of Pseudohynobius (Ps. flavomaculatus, Ps. guizhouensis, Ps. jinfo, Ps. kuankuoshuiensis and Ps. shuichengensis) based on x-ray computer tomography data for 18 specimens. Our results indicate that the cranium in each of these species has a combination of differences in morphology, proportions and articulation patterns in both dermal and endochondral bones. Our study establishes a range of intraspecific differences that will serve as organizing hypotheses for future studies as more extensive collections of these species become available. Morphological features in the cranium for terrestrial ecological adaptation in Hynobiidae are summarized. Based on the results, we also discuss the evolution and development of several potential synapomorphies of Hynobiidae, including features of the orbitosphenoid and articular.

Bid is involved in apoptosis induced by Chinese giant salamander iridovirus and contributes to the viral replication in an amphibian cell line.

Bid is a pro-apoptotic BH3-only member of the Bcl-2 superfamily that functions to link the extrinsic apoptotic pathway and the mitochondrial amplification loop of the intrinsic pathway. In this study, the expression and functions of Chinese giant salamander (Andrias davidianus) Bid (AdBid) were investigated. The AdBid cDNA sequence contains an open reading frame (ORF) of 576 nucleotides, encoding a putative protein of 191 aa. AdBid possesses the conserved BH3 interacting domain and shared 34-52% sequence identities with other amphibian Bid. mRNA expression of AdBid was most abundant in muscle. The expression level of AdBid in Chinese giant salamander muscle, kidney and spleen significantly increased after Chinese giant salamander iridovirus (GSIV) infection. Additionally, a plasmid expressing AdBid was constructed and transfected into the Chinese giant salamander muscle cell line (GSM cells). The morphology and cytopathic effect (CPE) and apoptotic process in AdBid over-expressed GSM cells was significantly enhanced during GSIV infection compared with that in control cells. Moreover, a higher level of the virus major capsid protein (MCP) gene copies and protein synthesis was confirmed in the AdBid over-expressed cells. These results indicated that AdBid played a positive role in GSIV induced apoptosis and the viral replication. This study may contribute to the better understanding on the infection mechanism of iridovirus-induced apoptosis.

Four immunoglobulin isotypes and IgD splice variants in urodele amphibians.

Until recently, different families of urodele amphibians were thought to express distinct subsets of immunoglobulin (Ig) isotypes. In this study, we explored cDNAs encoding Ig heavy-chains (H-chains) in three species of urodele amphibians. We found that Cynops pyrrhogaster, Pleurodeles waltl, and Ambystoma mexicanum each carry genes encoding four Ig H-chain isotypes, including IgM, IgY, IgD, and IgX, similar to those found in anuran amphibians. We also found that urodele IgDs have a long constant region similar to those found in anuran, reptiles, and bony fishes. We also found several putative IgD splice variants. Our findings indicated that P. waltl IgP is not a novel isotype but an IgD splice variant. Altogether, our findings indicate that IgD splice variants may be universally expressed among amphibian species.

Wide Distribution but Low Differentiation: Phylogenetic Relationships and Phylogeography of Hynobius nigrescens (Amphibia: Caudata).

We surveyed the genetic structure of Hynobius nigrescens Stejneger, 1907, a lentic breeding salamander widespread throughout montane and lowland regions of northeastern Japan. We performed a mitochondrial DNA analysis to explore intraspecific genetic variation and infer the evolutionary population history of H. nigrescens. Complete 1141 bp sequences of the mitochondrial cytochrome b gene were studied for 134 adult and larval individuals collected from 62 localities, encompassing the known range of the species. Hynobius nigrescens proved to be monophyletic, including two major clades (Clade II from southwestern Hokuriku and Clade I from all other localities). The latter clade comprises four well-supported and geographically structured subclades, which show genetic distances smaller than those seen in the widely sympatric species Hynobius lichenatus. Results of population statistical analyses indicated that Clade II of the westernmost range of H. nigrescens seems to have maintained a constant population size, while Clade I from most of the northeastern species range shows a tendency of recent population expansion, which is evident in Subclades I-A from the northernmost range and I-B from southern Tohoku to northern Kanto and eastern Chubu. In contrast, Subclades I-C from northeastern Chubu and Sado Is. and I-D from northwestern Chubu to Hokuriku seem to have been relatively stable in population size. Hynobius nigrescens differs greatly from other salamander species from northeastern Japan in its much more recent periods of genetic differentiation and its pattern thereof, and is suggested to be a young faunal element in this region.

Miniaturization, Genome Size, and Biological Size in a Diverse Clade of Salamanders.

AbstractGenome size (C-value) can affect organismal traits across levels of biological organization from tissue complexity to metabolism. Neotropical salamanders show wide variation in genome and body sizes, including several clades with miniature species. Because miniaturization imposes strong constraints on morphology and development and because genome size is strongly correlated with cell size, we hypothesize that body size has played an important role in the evolution of genome size in bolitoglossine salamanders. If this hypothesis is correct, then genome size and body size should be correlated in this group. Using Feulgen image analysis densitometry, we estimated genome sizes for 60 species of Neotropical salamanders. We also estimated the "biological size" of species by comparing genome size and physical body sizes in a phylogenetic context. We found a significant correlation between C-value and physical body size using optimal regression with an Ornstein-Uhlenbeck model and report the smallest salamander genome found to date. Our index of biological size showed that some salamanders with large physical body size have smaller biological body size than some miniature species and that several clades demonstrate patterns of increased or decreased biological size compared with their physical size. Our results suggest a causal relationship between physical body size and genome size and show the importance of considering the impact of both on the biological size of organisms. Indeed, biological size may be a more appropriate measure than physical size when considering phenotypic consequences of genome size evolution in many groups.

A new species of the genus Batrachuperus (Urodela: Hynobiidae) from Southwestern China.

A new stream salamander species, Batrachuperus daochengensis sp. nov., from southwestern China, is described herein based on morphological and molecular evidence. Molecular phylogeny derived from the mitochondrial gene together with previous nuclear data revealed that B. daochengensis sp. nov. is sister to B. yenyuanensis. The new species differs from all other species of the genus by the following combination of characters: brown horny epidermis on tips of fingers and toes absent; tubercles on palms and soles absent; costal grooves 12; dorsal brown, mottled with blackish spots; fingers 2-3-4-1 in order of decreasing length; tips of longest digits of fore- and hindlimbs largely separated by one to two costal spaces when adpressed towards each other along sides of body. The new species is currently known in the central and southern Shaluli Mountains in southwestern China.

A Triassic stem-salamander from Kyrgyzstan and the origin of salamanders.

The origin of extant amphibians remains largely obscure, with only a few early Mesozoic stem taxa known, as opposed to a much better fossil record from the mid-Jurassic on. In recent time, anurans have been traced back to Early Triassic forms and caecilians have been traced back to the Late Jurassic Eocaecilia, both of which exemplify the stepwise acquisition of apomorphies. Yet the most ancient stem-salamanders, known from mid-Jurassic rocks, shed little light on the origin of the clade. The gap between salamanders and other lissamphibians, as well as Paleozoic tetrapods, remains considerable. Here we report a new specimen of Triassurus sixtelae, a hitherto enigmatic tetrapod from the Middle/Late Triassic of Kyrgyzstan, which we identify as the geologically oldest stem-group salamander. This sheds light not only on the early evolution of the salamander body plan, but also on the origin of the group as a whole. The new, second specimen is derived from the same beds as the holotype, the Madygen Formation of southwestern Kyrgyzstan. It reveals a range of salamander characters in this taxon, pushing back the rock record of urodeles by at least 60 to 74 Ma (Carnian-Bathonian). In addition, this stem-salamander shares plesiomorphic characters with temnospondyls, especially branchiosaurids and amphibamiforms.

Vastly underestimated species richness of Amazonian salamanders (Plethodontidae: Bolitoglossa) and implications about plethodontid diversification.

We present data showing that the number of salamander species in Amazonia is vastly underestimated. We used DNA sequences of up to five genes (3 mitochondrial and 2 nuclear) of 366 specimens, 189 corresponding to 89 non-Amazonian nominal species and 177 Amazonian specimens, including types or topotypes, of eight of the nine recognized species in the region. By including representatives of all known species of Amazonian Bolitoglossa, except for one, and 73% of the currently 132 recognized species of the genus, our dataset represents the broadest sample of Bolitoglossa species, specimens, and geographic localities studied to date. We performed phylogenetic analyses using parsimony with tree-alignment and maximum likelihood (ML) with similarity alignment, with indels as binary characters. Our optimal topologies were used to delimit lineages that we assigned to nominal species and candidate new species following criteria that maximize the consilience of the current species taxonomy, monophyly, gaps in branch lengths, genetic distances, and geographic distribution. We contrasted the results of our species-delimitation protocol with those of Automated Barcode Gap Discovery (ABGD) and multi-rate Poisson Tree Processes (mPTP). Finally, we inferred the historical biogeography of South American salamanders by dating the trees and using dispersal-vicariance analysis (DIVA). Our results revealed a clade including almost all Amazonian salamanders, with a topology incompatible with just the currently recognized nine species. Following our species-delimitation criteria, we identified 44 putative species in Amazonia. Both ABGD and mPTP inferred more species than currently recognized, but their numbers (23-49) and limits vary. Our biogeographic analysis suggested a stepping-stone colonization of the Amazonian lowlands from Central America through the Chocó and the Andes, with several late dispersals from Amazonia back into the Andes. These biogeographic events are temporally concordant with an early land bridge between Central and South America (~10-15 MYA) and major landscape changes in Amazonia during the late Miocene and Pliocene, such as the drainage of the Pebas system, the establishment of the Amazon River, and the major orogeny of the northern Andes.

A new small-sized stem salamander from the Middle Jurassic of Western Siberia, Russia.

Salamanders (Caudata) are one of the three modern groups of amphibians known from the Middle Jurassic. The early stages of evolution of these amphibians are still poorly known, especially for stem taxa of Jurassic age. A new small-sized stem salamander, Egoria malashichevi gen. et sp. nov., from the Middle Jurassic (Bathonian) Itat Formation of the Berezovsk Quarry locality in Western Siberia, Russia, is described on the basis of isolated vertebrae, including an atlas centrum and a fragmentary trunk vertebra centrum previously referred to an undescribed salamander taxon ("Berezovsk salamander A"). The new taxon is diagnosed by the following unique combination of vertebral characters: atlantal anterior cotyles with elliptical anterior outline, located at an angle of approximately 135-137 degrees to each other; wide posterior portion of the atlantal centrum; ossified portion of the intercotylar tubercle represented by dorsal and ventral lips; absence of a deep depression on the ventral surface of the atlantal centrum; absence of pronounced ventrolateral ridges on the atlas; absence of spinal nerve foramina; presence of a pitted texture on the ventral and lateral surfaces of the centra and lateral surfaces neural arch pedicels; presence of a short atlantal neural arch with its anterior border situated behind the level of the anterior cotyles; short trunk vertebrae; and upper transverse process (= diapophysis) larger than lower transverse process (= parapophysis) on the trunk vertebrae; notochordal canal opens in the upper half of the cotyle (= the lower portion of the centrum is more massive and less compact than the upper portion). The microanatomical organization of the atlas and trunk vertebrae is characterized by the presence of inner cancellous endochondral bone. The small body size (about 180-215 mm) of Egoria malashichevi gen. et sp. nov. indicates that that not all stem salamanders were large neotenic forms (up to 550-600 mm in Urupia and Marmorerpeton) and hints at a broader ecological role for stem salamanders.

The last Palaeoproteus (Urodela: Batrachosauroididae) of Europe.

The Batrachosauroididae are an enigmatic group of salamanders known from the Cretaceous and Tertiary of North America and Europe. In Europe, the family is known only by two species of the genus Palaeoproteus. The genus has limited distribution in Western and Central Europe. In the present paper, we describe a new species, Palaeoproteus miocenicus, from the early late Miocene (11-9 Ma) of Austria and Ukraine, representing the youngest record of the family Batrachosauroididae from the Neogene of Europe. The new species differs from the Paleogene representatives of the genus by 12 characters, including large body size, the long anterior extension of the Meckelian groove and the size and shape of the odontoid process on the dentary. The µCT scanning of bones of the new species revealed novel features (e.g. anterior extension of Meckelian groove, interconnected network of canals and small cavities in atlas) observable only in this species. P. miocenicus inhabited aquatic environments, which existed under wet climatic conditions with mean annual precipitation higher than 900 mm. The new species expands the temporal range of the genus by at least 30 million years and enlarges the palaeogeographic distribution of the genus into Eastern Europe.

Salamanders and other amphibians are aglow with biofluorescence.

Biofluorescence is the absorption of electromagnetic radiation (light) at one wavelength followed by its reemission at a lower energy and longer wavelength by a living organism. Previous studies have documented the widespread presence of biofluorescence in some animals, including cnidarians, arthropods, and cartilaginous and ray-finned fishes. Many studies on biofluorescence have focused on marine animals (cnidarians, cartilaginous and ray-finned fishes) but we know comparatively little about the presence of biofluorescence in tetrapods. We show for the first time that biofluorescence is widespread across Amphibia, with a focus on salamanders (Caudata), which are a diverse group with a primarily Holarctic distribution. We find that biofluorescence is not restricted to any particular family of salamanders, there is striking variation in their fluorescent patterning, and the primary wavelengths emitted in response to blue excitation light are within the spectrum of green light. Widespread biofluorescence across the amphibian radiation is a previously undocumented phenomenon that could have significant ramifications for the ecology and evolution of these diverse and declining vertebrates. Our results provide a roadmap for future studies on the characterization of molecular mechanisms of biofluorescence in amphibians, as well as directions for investigations into the potential impact of biofluorescence on the visual ecology and behavior of biofluorescent amphibians.

Detecting Mosaic Patterns in Macroevolutionary Disparity.

Evolutionary biologists have long sought to understand the full complexity in pattern and process that shapes organismal diversity. Although phylogenetic comparative methods are often used to reconstruct complex evolutionary dynamics, they are typically limited to a single phenotypic trait. Extensions that accommodate multiple traits lack the ability to partition multidimensional data sets into a set of mosaic suites of evolutionarily linked characters. I introduce a comparative framework that identifies heterogeneity in evolutionary patterns across large data sets of continuous traits. Using a model of continuous trait evolution based on the differential accumulation of disparity across lineages in a phylogeny, the approach algorithmically partitions traits into a set of character suites that best explains the data, where each suite displays a distinct pattern in phylogenetic morphological disparity. When applied to empirical data, the approach revealed a mosaic pattern predicted by developmental biology. The evolutionary distinctiveness of individual suites can be investigated in more detail either by fitting conventional comparative models or by directly studying the phylogenetic patterns in disparity recovered during the analysis. This framework can supplement existing comparative approaches by inferring the complex, integrated patterns that shape evolution across the body plan from disparate developmental, morphometric, and environmental sources of phenotypic data.

Phylogenomic data reveal reticulation and incongruence among mitochondrial candidate species in Dusky Salamanders (Desmognathus).

Gene flow between evolutionarily distinct lineages is increasingly recognized as a common occurrence. Such processes distort our ability to diagnose and delimit species, as well as confound attempts to estimate phylogenetic relationships. A conspicuous example is Dusky Salamanders (Desmognathus), a common model-system for ecology, evolution, and behavior. Only 22 species are described, 7 in the last 40 years. However, mitochondrial datasets indicate the presence of up to 45 "candidate species" and multiple paraphyletic taxa presenting a complex history of reticulation. Some authors have even suggested that the search for species boundaries in the group may be in vain. Here, we analyze nuclear and mitochondrial data containing 161 individuals from at least 49 distinct evolutionary lineages that we treat as candidate species. Concatenated and species-tree methods do not estimate fully resolved relationships among these taxa. Comparing topologies and applying methods for estimating phylogenetic networks, we find strong support for numerous instances of hybridization throughout the history of the group. We suggest that these processes may be more common than previously thought across the phylogeography-phylogenetics continuum, and that while the search for species boundaries in Desmognathus may not be in vain, it will be complicated by factors such as crypsis, parallelism, and gene-flow.